A Theoretical Perspective on the Photochemistry of Boron-Nitrogen Lewis Adducts.

Boron-Nitrogen (B-N) Lewis adducts form a versatile family of compounds with numerous applications in functional molecules. Despite the growing interest in this family of compounds for optoelectronic applications, little is currently known about their photophysics and photochemistry. Even the electronic absorption spectrum of ammonia borane, the textbook example of a B-N Lewis adduct, is unavailable. Given the versatility of the light-induced processes exhibited by these molecules, we propose in this work a detailed theoretical study of the photochemistry and photophysics of simple B-N Lewis adducts. We used advanced techniques in computational photochemistry to identify and characterize the possible photochemical pathways followed by ammonia borane and extended this knowledge to the substituted B-N Lewis adducts pyridine-borane and pyridine-boric acid. The photochemistry observed for this series of molecules allows us to extract qualitative rules to rationalize the light-induced behavior of more complex B-N-containing molecules.

On the fluorescence enhancement of arch neuronal optogenetic reporters.

The lack of a theory capable of connecting the amino acid sequence of a light-absorbing protein with its fluorescence brightness is hampering the development of tools for understanding neuronal communications. Here we demonstrate that a theory can be established by constructing quantum chemical models of a set of Archaerhodopsin reporters in their electronically excited state. We found that the experimentally observed increase in fluorescence quantum yield is proportional to the computed decrease in energy difference between the fluorescent state and a nearby photoisomerization channel leading to an exotic diradical of the protein chromophore. This finding will ultimately support the development of technologies for searching novel fluorescent rhodopsin variants and unveil electrostatic changes that make light emission brighter and brighter.

A Theoretical Perspective on the Actinic Photochemistry of 2-Hydroperoxypropanal.

The photochemical reactions triggered by the sunlight absorption of transient volatile organic compounds in the troposphere are notoriously difficult to characterize experimentally due to the unstable and short-lived nature of these organic molecules. Some members of this family of compounds are likely to exhibit a rich photochemistry given the diversity of functional groups they can bear. Even more interesting is the photochemical fate of volatile organic compounds bearing more than one functional group that can absorb light─this is the case, for example, of α-hydroperoxycarbonyls, which are formed during the oxidation of isoprene. Experimental observables characterizing the photochemistry of these molecules like photoabsorption cross-sections or photolysis quantum yields are currently missing, and we propose here to leverage a recently developed computational protocol to predict in silico the photochemical fate of 2-hydroperoxypropanal (2-HPP) in the actinic region. We combine different levels of electronic structure methods─SCS-ADC(2) and XMS-CASPT2─with the nuclear ensemble approach and trajectory surface hopping to understand the mechanistic details of the possible nonradiative processes of 2-HPP. In particular, we predict the photoabsorption cross-section and the wavelength-dependent quantum yields for the observed photolytic pathways and combine them to determine in silico photolysis rate constants. The limitations of our protocol and possible future improvements are discussed.

Quantum dynamics with curvilinear coordinates: models and kinetic energy operator.

In order to simplify the numerical solution of the time-dependent or time-independent Schrödinger equations associated with atomic and molecular motions, the use of well-adapted coordinates is essential. Usually, this set of curvilinear coordinates leads to a Hamiltonian operator that is as separable as possible. Although their corresponding kinetic energy operator (KEO) expressions can be derived analytically for small systems or special kinds of coordinates, a numerical and exact approach allows one to compute them in terms of sophisticated curvilinear coordinates. Furthermore, the numerical approach enables one to easily define reduced-dimensionality or constrained models. We present here a recent implementation of this numerical approach that allows nested coordinate transformations, therefore leading to great flexibility in the definition of the curvilinear coordinates. Furthermore, this implementation has no limitations in terms of numbers of atoms or coordinate transformations. The quantum dynamics of the cis-trans photoisomerization of part of the retinal chromophore illustrates the construction of the coordinates and KEO part of a three-dimensional model. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.

Calculating Photoabsorption Cross-Sections for Atmospheric Volatile Organic Compounds.

Characterizing the photochemical reactivity of transient volatile organic compounds (VOCs) in our atmosphere begins with a proper understanding of their abilities to absorb sunlight. Unfortunately, the photoabsorption cross-sections for a large number of transient VOCs remain unavailable experimentally due to their short lifetime or high reactivity. While structure-activity relationships (SARs) have been successfully employed to estimate the unknown photoabsorption cross-sections of VOCs, computational photochemistry offers another promising strategy to predict not only the vertical electronic transitions of a given molecule but also the width and shape of the bands forming its absorption spectrum. In this work, we focus on the use of the nuclear ensemble approach (NEA) to determine the photoabsorption cross-section of four exemplary VOCs, namely, acrolein, methylhydroperoxide, 2-hydroperoxy-propanal, and (microsolvated) pyruvic acid. More specifically, we analyze the influence that different strategies for sampling the ground-state nuclear density-Wigner sampling and ab initio molecular dynamics with a quantum thermostat-can have on the simulated absorption spectra. We highlight the potential shortcomings of using uncoupled harmonic modes within Wigner sampling of nuclear density to describe flexible or microsolvated VOCs and some limitations of SARs for multichromophoric VOCs. Our results suggest that the NEA could constitute a powerful tool for the atmospheric community to predict the photoabsorption cross-section for transient VOCs.

Caveat when using ADC(2) for studying the photochemistry of carbonyl-containing molecules.

Several electronic-structure methods are available to study the photochemistry and photophysics of organic molecules. Among them, ADC(2) stands as a sweet spot between computational efficiency and accuracy. As a result, ADC(2) has recently seen its number of applications booming, in particular to unravel the deactivation pathways and photodynamics of organic molecules. Despite this growing success, we demonstrate here that care has to be taken when studying the nonradiative pathways of carbonyl-containing molecules, as ADC(2) appears to suffer from a systematic flaw.

Relaxation dynamics through a conical intersection: Quantum and quantum-classical studies.

We study the relaxation process through a conical intersection of a photo-excited retinal chromophore model. The analysis is based on a two-electronic-state two-dimensional Hamiltonian developed by Hahn and Stock [J. Phys. Chem. B 104 1146 (2000)] to reproduce, with a minimal model, the main features of the 11-cis to all-trans isomerization of the retinal of rhodopsin. In particular, we focus on the performance of various trajectory-based schemes to nonadiabatic dynamics, and we compare quantum-classical results to the numerically exact quantum vibronic wavepacket dynamics. The purpose of this work is to investigate, by analyzing electronic and nuclear observables, how the sampling of initial conditions for the trajectories affects the subsequent dynamics.

Quantum and Quantum-Classical Studies of the Photoisomerization of a Retinal Chromophore Model.

We report an in-depth analysis of the photo-induced isomerization of the 2-cis-penta-2,4-dieniminium cation: a minimal model of the 11-cis retinal protonated Schiff base chromophore of the dim-light photoreceptor rhodopsin. Based on recently developed three-dimensional potentials parametrized on ab initio multi-state multi-configurational second-order perturbation theory data, we perform quantum-dynamical studies. In addition, simulations based on various quantum-classical methods, among which Tully surface hopping and the coupled-trajectory approach derived from the exact factorization, allow us to validate their performance against vibronic wavepacket propagation and, therefore, a purely quantum treatment. Quantum-dynamics results uncover qualitative differences with respect to the two-dimensional Hahn-Stock potentials, widely used as model potentials for the isomerization of the same chromophore, due to the increased dimensionality and three-mode correlation. Quantum-classical simulations show, instead, that three-dimensional model potentials are capable of capturing a number of features revealed by atomistic simulations and experimental observations. In particular, a recently reported vibrational phase relationship between double-bond torsion and hydrogen-out-of-plane modes critical for rhodopsin isomerization efficiency is correctly reproduced.

On the Theoretical Determination of Photolysis Properties for Atmospheric Volatile Organic Compounds.

Volatile organic compounds (VOCs) are ubiquitous atmospheric molecules that generate a complex network of chemical reactions in the troposphere, often triggered by the absorption of sunlight. Understanding the VOC composition of the atmosphere relies on our ability to characterize all of their possible reaction pathways. When considering reactions of (transient) VOCs with sunlight, the availability of photolysis rate constants, utilized in general atmospheric models, is often out of experimental reach due to the unstable nature of these molecules. Here, we show how recent advances in computational photochemistry allow us to calculate in silico the different ingredients of a photolysis rate constant, namely, the photoabsorption cross-section and wavelength-dependent quantum yields. The rich photochemistry of tert-butyl hydroperoxide, for which experimental data are available, is employed to test our protocol and highlight the strengths and weaknesses of different levels of electronic structure and nonadiabatic molecular dynamics to study the photochemistry of (transient) VOCs.

Two-State, Three-Mode Parametrization of the Force Field of a Retinal Chromophore Model.

In recent years, the potential energy surfaces of the penta-2,4-dieniminium cation have been investigated using several electronic structure methods. The resulting pool of geometrical, electronic, and energy data provides a suitable basis for the construction of a topographically correct analytical model of the molecule force field and, therefore, for a better understanding of this class of molecules, which includes the chromophore of visual pigments. In the present contribution, we report the construction of such a model for regions of the force field that drive the photochemical and thermal isomerization of the central double bound of the cation. While previous models included only two modes, it is here shown that the proposed three-mode model and corresponding set of parameters are able to reproduce the complex topographical and electronic structure features seen in electronically correlated data obtained at the XMCQDPT2//CASSCF/6-31G* level of theory.